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Based on the wormlike chain model, a coarse-grained description of the nonlinear dynamics of 
a weakly bending semiflexible polymer is developed. By means of a multiple scale perturbation 
analysis, a length-scale separation inherent to the weakly-bending limit is exploited to reveal the 
deterministic nature of the spatio-temporal relaxation of the backbone tension and to deduce the 
corresponding coarse-grained equation of motion. From this partial integro-differential equation, 
some detailed analytical predictions for the non-linear response of a weakly bending polymer are 
derived in an accompanying paper (Part II [13). 
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I. INTRODUCTION 



The laws of Brownian motion have played the role of 
a mediator between the apparently smooth determin- 
istic dynamics on a macroscopic scale and the micro- 
scopic molecular chaos since their discovery a century 
ago . They are pivotal to our understanding of a broad 
class of animate and inanimate soft condensed matter 
systems that owe their characteristic softness to low- 
dimensional and strongly fluctuating meso-scale struc- 
tures such as polymeric networks and membranes 0- 
Conversely, these systems are well suited to study how 
complex deterministic dynamics on a macro- or meso- 
scale emerges from the underlying stochastic differential 
equations 

Take, for example, a stiff polymer like actin that is 
suddenly stretched by strong forces applied at its ends. 
Or, conversely, consider a polymer that is held in a vir- 
tually straight conformation and suppose that the forces 
at its ends are suddenly released. These two paradig- 
matic experimental setups, which we call Pulling and 
Release, are illustrated in Fig. ^ How will the end-to- 
end distance of the polymer relax to its new equilibrium 
value? Given the manifestly stochastic underlying dy- 
namics, which for Release is exclusively driven by thermal 
forces, it is not immediately obvious that the initial con- 
traction or stretching dynamics should obey a determin- 
istic law, as tacitly assumed by several heuristic deriva- 
tions Indeed, these studies, which predicted 
a variety of interesting new dynamic scaling regimes, ar- 
rived at partially contradicting results 0,0. Despite con- 
siderable experimental, theoretical and numerical work, 
both relaxation laws (for Pulling and Release) remained 
controversial for some time. 
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It therefore appears worthwhile to undertake a detailed 
mathematical derivation of the meso-scale dynamic equa- 
tions that govern the nonlinear dynamics of semiflexible 
polymers "from first principles" Q; i.e. from the under- 
lying stochastic differential equations of motion. In a re- 
cent Letter, we have outlined such a systematic approach 
that resolves the aforementioned theoretical problems, 
together with some of its consequences The present 
contribution offers a more comprehensive discussion. In 
Part I, an effective coarse-grained meso-scale description 
of the dynamics of a semiflexible polymer is developed by 
means of a multiple scale theory from the stochastic dif- 
ferential equations of motion. Our detailed analysis also 
reveals the limits of validity of the deterministic meso- 
scopic description and shows how to deal with subtle end 
effects that may in some cases mask the non-trivial pre- 
dictions for certain observables. Building on this general 
framework, the theory is elaborated for the specific prob- 
lems of Pulling and Release in Part II |lCll |. which pro- 
vides a template for the future analysis of a variety of re- 
lated problems with somewhat different boundary/initial 
conditions [Ullllll. Thereby, we corroborate the im- 
portance of a regime of homogeneous tension relaxation, 
which generally occurs in i?e/ease-experiments, and es- 
tablish the complete crossover scenario between the vari- 
ous intermediate scaling regimes. Additionally, in Part II 
some consequences for common observables are worked 
out in detail to facilitate experimental verification of the 
theory. 

Before entering a detailed quantitative analysis, it 
seems useful to summarize the main ideas on a quali- 
tative level in order to make the remainder more easily 
accessible. A characteristic property of semiflexible poly- 
mers and many other fluctuating meso-scale structures in 
soft condensed matter, is their reduced dimensionality or 
slender shape. It entails the presence of thermally ex- 
cited transverse fluctuations of an essentially inextensi- 
ble backbone. Returning to the above example, an actin 
filament is much more susceptible to bending undula- 
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FIG. 1: Two basic examples of dynamic force-extension ex- 
periments. Pulling: a weakly bending polymer in equilibrium 
is suddenly pulled longitudinally by two external forces f . 
Release: a pre-stretched polymer is suddenly released. 



tions than to stretching or compressing its backbone. An 
analogous statement holds holds for other biopolymers or 
two-dimensional locally flat objects like membranes and 
surfaces. This suggests to idealize these structures as un- 
dulating inextensible manifolds |l4j . To be specific, we 
focus for the following on the case of a single semiflexi- 
ble polymer in solution, which seems to be the simplest 
paradigmatic example of the more general soft-matter 
meso-structures alluded to above. Its equilibrium me- 
chanical properties and conformational statistics have by 
now been thoroughly studied theoretically and experi- 
mentally [n m. Both are well understood in terms 
of the sclf-affinc roughness of the equilibrium contour 
fluctuations within the so-called wormlike chain model, 
which idealizes the polymer as an inextensible space 
curve with an energetic cost for bending [T^ . Here, we 
are primarily interested in the transient non-equilibrium 
stretching and contraction dynamics of such a wormlike 
chain, i.e., in the much less studied problem of how a 
scmiflexiblc polymer relaxes to equilibrium after a sud- 
den drastic change in its boundary conditions. This ques- 
tion is of considerable fundamental and practical interest 
alike, e.g. for single-molecule manipulations [T^ ITfll l2nl | 
and for understanding and controlling the dynamic re- 
sponse of polymer solutions and networks such as those 
constituting the cytoskeleton of biological cells [U H^l . 

Because of the inextensible backbone, the stretching 
or contraction dynamics of a wormlike chain is entirely 
due to a spatio-temporal re-organization of the contour 
length stored in the transverse thermal wrinkles. It is 
governed by the dynamic backbone tension /(s,t), which 
is the force that holds the backbone together. The dy- 
namics is always assumed to be strongly overdamped by 
solvent friction, which can be decomposed into transverse 
and longitudinal components with respect to the local 
tangent, in view of the locally rod-like conformation of 
the polymer. At flrst sight, one might suppose that for 
small enough undulations one can resort to a formulation 
of the dynamics solely in terms of transverse modes for 
which only transverse friction matters. Actually, for the 
equilibrium dynamics in the absence of external forces, 
the conclusions based on such an assumption are in ac- 



cord with a series of experimental data |23ll24ll25ll26j . At 

a second thought, considering the constraint of inexten- 
sibility, it is far from obvious how the interplay between 
transverse and longitudinal friction limits the relaxation 
after a sudden application or release of external forces 
as e.g. in Release and Pulling. It is the major objective 
of the present study to resolve this puzzle for the case of 
longitudinal forces, while the somewhat more complex is- 
sue of the nonlinear transverse response will be addressed 
in a separate contribution • 

As a cornerstone of our derivation, we establish a dy- 
namic scale separation between the scales where trans- 
verse and longitudinal friction reign, namely the trans- 
verse and longitudinal dynamic correlation lengths i±{t) 
and £\\{t), respectively. We demonstrate that £±{t) <C 
£\\ (t) holds at any time in the limit of a weakly bending 
rod. As a central result, we obtain that the tension varies 
to leading order only on the large scale £\\{t), over which 
the short-wavelength transverse undulations that domi- 
nate the dissipation on the scale £±{t) arc self-averaging. 
The short-scale transverse fluctuations may thus be said 
to provide an effective local backbone compliance for the 
large-scale longitudinal dynamics. In this way, transverse 
and longitudinal dynamics are effectively decoupled, and 
the microscopic stochastic differential equations can be 
reduced, in a controlled way, to a deterministic (integro- 
differential) equation for the coarse-grained meso-scale 
dynamics. 

With some care, most of the predictions that emanate 
from this reduced description can qualitatively be ob- 
tained from a relatively simple scaling analysis, to which 
we will occasionally resort in order to promote the in- 
tuitive physical understanding of our analysis. In fact, 
not only physical insights but also several interesting for- 
mal predictions were originally obtained on the basis of 
simple scaling arguments. However, due to some sub- 
tleties and pitfalls, the results available in the literature 
remained somewhat contradictory. The following devel- 
opment, in particular Part I of this article scries, aims at 
settling the corresponding issues conclusively by employ- 
ing a controlled perturbation theory instead of evoking 
scaling assumptions. 

As suggested by the introductory examples of dynamic 
force-extension experiments, exemplary realizations of 
the meso-scale dynamics occur in response to highly lo- 
calized forces, which we generally (though not invari- 
ably) assume to be applied abruptly at the boundaries. 
They initiate universal self-similar relaxation processes 
that spread through the polymer, resulting in character- 
istic power-law signatures, so-called intermediate asymp- 
totics j23l, in various observables. These will be derived 
analytically in Part II of this article, which is devoted 
to solving the effective deterministic meso-scale equations 
for the tension, obtained below. There, we will also con- 
sider the consequences of the tension dynamics on perti- 
nent observables like the (projected) end-to-end distance 
and the novel experimental perspectives brought up b" 
our analysis. For the impatient reader, our Letter 
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may serve as a quick guide to our main arguments and 
results. 

The outline of the present part (Part I), is as fol- 
lows. We develop the systematic coarse grained descrip- 
tion of tension dynamics in stiff, respectively, stretched 
semiflexible polymers that emerges from an appropri- 
ate small gradient expansion of the dynamical worm- 
like chain model (Sec. An ordinary perturbation 
expansion leading to linear first order equations of mo- 
tion (Sec. Illl|l turns out to be restricted to times when 
the tension has already relaxed to its equilibrium value 
fSec. IIV|) . The actual tension dynamics on shorter times 
is resolved by a stochastic multiple scale perturbation 
theory (Sec. 0, which is based on a dynamical length 
scale separation. As a major result, we obtain a rigorous 
deterministic partial integro-differcntial equation (PIDE) 
that describes the long wave-length (-all time-) dynamics 
of the tension. 



II. DYNAMICAL WORMLIKE CHAIN MODEL 



sibility constraint, Eq. (QJ, has to be enforced as a rigid 
constraint. 

The elastic force (per arclength) g^x derives from 
Eq. I^J by a functional derivative, 



5oi(s,i) = - 



5n 



WLC 



5r 



(3) 



As indicated, the contour variations 5r used to determine 
the functional derivative on the right-hand-side have to 
respect the local inextensibility constraint, Eq. (QJ. The 
variational calculation, detailed in App. ^ yields 



r' X (Kr"' + /ei) = 0, 



(4) 



where /d (s) = J^dsgci{s) is the spatial integral over the 
elastic force density defined in Eq. 

According to the implicit equation for fd , Eq. Q , the 
force K r'" + fd has a vanishing component transverse to 
the local tangent r' . Equivalently, we may require that 
both vectors are proportional to each other, 



At low Reynolds numbers, the dynamics of a polymer 
is determined by the balance of elastic forces, friction, 
and stochastic forces. We will briefly motivate how these 
forces are modeled in the usual stochastic description of 
the Brownian motion of a semiflexible polymer, leading 
to Eq. below, which is the basis of our subsequent 
analysis. 

The natural model for the description of the elas- 
tic properties of a semiflexible polymer arises from the 
idea to regard the polymer as a thin cylinder consist- 
ing of a homogeneous elastic material. In the slender 
limit, where the ratio of thickness over total length ap- 
proaches zero, the deformation modes of the cylinder 
become much stiffer in the axial direction ("phonons") 
than in directions transverse to the cylinder axis ( "bend- 
ing modes") Consequently, such a slender rod or 
"thread" subject to external forces merely allows for 
bending deformations, while its contour length is to a 
good approximation locally conserved. 

These features are idealized in the wormlike chain 
model, where the polymer is at any time t represented as 
an inextensible space curve r{s,t) parameterized by the 
arc length s, i.e. the tangents have to obey the constraint 

r'2 = 1 . (1) 

Here, we have introduced the shorthand notation r' = 
dr{s, t)/ds. The efi^ective free energy Hwlc of a particu- 
lar conformation is only due to bending (curvature) and 
is given by 

■HwLC^^ iysr"\ (2) 



^^"'+/ci =/r', (5) 

where the proportionality factor f{s) has dimensions of 
a force and, in general, depends on the arc length. A 
spatial derivative then yields a direct expression for the 
elastic force density, 

gd = -^r"" + {fr'y . (6) 

The apparent simplification with respect to Eq. I^J is 
somewhat deceiving, since Eq. ^ still contains the un- 
known function /(s), which has to be fixed by the local 
arc-length constraint, Eq. QJ. The new force field /(s) 
has however a very direct and intuitive physical inter- 
pretation. It is the local line tension |29l | that "ties the 
polymer together" and thereby enforces the inextensibil- 
ity condition. 

Much of the following deals with the dynamics of this 
local line tension, which turns out to be closely related to 
the dynamics of the local excess length stored in undu- 
lations, defined in Eq. IjlSfl below, to which we will refer 
to as stored length. Time-dependence is introduced into 
the description by requiring the clastic force density gd 
to be balanced by the dynamic friction (per arc length) 
with the solvent, 

gfr{s,t) ^ -<^{r,t)dtr{s,t) (free draining) . (7) 

At any arclength s and time t, the friction matrix ^ with 
elements Qj can be decomposed into its transverse and 
longitudinal components with respect to the local tangent 
r'{s,t), 

= (1 - r' ® r') + C|| r' (g) r'] . (8) 



where k is the bending stiffness. To assure that the in- 
tegrand is the square of the local curvature, the inexten- 



The constants and can to a first approximation 
be estimated by the friction coefficients (per length) for 
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transverse respectively longitudinal motion of a rigid 
slender rod in a solvent of viscosity r] |30l |. 



a = 2C|| 47rr/ 



(9) 



A more sophisticated analysis would consider logarithmic 
corrections Is^ l33l | to account for the dynamic cou- 
pling of distant chain segments s and s' via long-ranged 
hydrodynamic interactions [s^, 01 ■ While such logarith- 
mic factors may sometimes be crucial in a quantitative 
comparison with some experiments [23, , their (fea- 
sible) implementation is not of primary interest to our 
present discussion, so that we chose to dismiss them for 
greater clarity of the presentation. 

On the same level of approximation, the force balance 
of elastic and frictional forces can be extended by adding 
thermal white noise fully characterized in terms of its 
mean and variance. 



0, 



(10a) 

s')S{t-t') (10b) 



with the angular brackets indicating an ensemble aver- 
age. The strength of the noise correlations in Eq. IllObll i s 
dictated by the Fluctuation-Dissipation theorem 
which assures that the steady state of the stochastic equa- 
tions of motion correspond to thermal equilibrium (for 
complications in simulations of discrete bead-rod chains, 
see Ref. [aSSIli)- 

Upon using Eqs. 1001, the balance of elastic, friction 
and stochastic forces, = gc\+ Sir + ^, takes the form of 
an equation of motion, 



C9tr = -Atr-"" + (/r')'+^ . 



(11) 



The partial differential equation Eq. Hll() . the arc length 
constraint Eq. (0) and the Gaussian noise-correlation 
comprise a complete stochastic description of the 
Brownian dynamics in the free draining limit. 



III. LINEARIZED STOCHASTIC DYNAMICS 

The nonlinear stochastic dynamics of a wormlike chain, 
represented by Eq. Hll|l . in combination with the inexten- 
sibility constraint, Eq. J^l, is hard to analyze in general. 
The difficulties are largely due to the calculation of the 
line tension /(s,t), which has to enforce the local inex- 
tensibility constraint, Eq. (0). There have been attempts 
to relax the local constraint and merely enforce local in- 
extensibility on average ji^, 0, E^l • This corresponds to 
replacing the field /(s,t) by a (spatially averaged) mean 
field. However those approaches fail to describe semiflex- 
ible polymers on local scales, where tension fluctuations 
(in time and space) are important 

Our analytical approach is instead based on the weakly 
bending limit, in which the polymer conformation is ap- 
propriately described by its deviations from a straight 




FIG. 2: Typical conformations of a freely fluctuating stiff 
polymer with suppressed rotations. In order to generate these 
conformations, which serve illustrative purposes only, we have 
represented each conformation as a linear superposition of the 
first 10 modes that solve the linearized equation of motion. 
For each realization, the mode amplitudes were then drawn 
randomly according to their Gaussian distribution in equilib- 



line. Accordingly, we parameterize the polymer's space 
curve by 



1 



(cf. Fig. 13, 



where r±{s) and r||(s) are the transverse and longitudi- 
nal displacements at arc length s. Such parameterization 
is only useful if the gradients of the transverse displace- 
ments are small everywhere along the contour, 



r'l = 0{e) < 1 



s e (0,L) . 



(12) 



The weakly bending limit, as defined in Eq. (|12|l . assumes 
that there is a small parameter e that controls the poly- 
mer roughness uniformly along the contour. In the case 
of a stiff polymer, on which we focus if otherwise stated, 
the small parameter is provided by the ratio 



LkeT 



< 1 



(stiff polymer) (13) 



of length L over persistence length ip. Alternatively, a 
weakly bending conformation may be realized by apply- 
ing a static external stretching force f larger than the 
internal characteristic force scale ksT/lp. In this case, 
r'l ^ y^kBT/ipf liil and we may identify 



ksT/epf (stretched polymer). (14) 



Eq. (fT^ allows us to expand the dynamical equations 
of motion in terms of small gradients. We start with 
the inextensibility constraint, which "slaves" the higher 
order longitudinal displacements to the transverse ones. 
After resolving Eq. ^ for r'^^ and expanding the square 
root, the local constraint takes the simple form 



(15) 



This entails that the parameter e is a measure for the 
reduction of the longitudinal extension 



=L-r||(L)+r||(0) 



(16) 



of the polymer due to the presence of thermal undula- 
tions: An arclength integral of Eq. (|15|) shows that is 
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smaller then the total contour length L by an amount of 
the order 

L-i?|| =r||(L)-r||(0) = O(eL). (17) 

We may think of the length difference in Eq. (|17() as being 
"stored" in undulations. The distribution of this excess 
length along the filament is, according to Eq. 1)15(1 . de- 
scribed by the function 

eM^\r'l{s4) = 0{e), (18) 

which will have central importance in our analysis. It is 
the fraction of the contour length that is at arclength s 
and time t locally stored in undulations. 

After taking a "spatial" (i.e. arc-length) derivative of 
its longitudinal part, we expand the equation of motion, 
Eq. ijTTI), to order 0{r'l) = 0{e) and obtain 

dtv^ = -rr + (/rl)'+fl±+€^ (19a) 
Qdtr\ = (C-l)(rl9,r^)'-rp-r + (/rf|)" 

-flil-^ll- (19b) 

Here, we have neglected terms of order 0{f^^'^) and made 
the following choice of units: time and tension, respec- 
tively, are rescaled according to 

t C±t/K (20) 
/ '^Z • (21) 

This corresponds to setting k = (± = 1 and ( = 1/2 = C|| ■ 
As a consequence all variables represent powers of length, 
e.g., t and / are a length''^ and a length"^, respectively. 
In Eq. H19() we have further allowed for an external |5J| 
force density g — (fl±,0||), which is a length"'^, that 
may for instance represent the effect of a solvent flow, an 
optical/magnetical tweezer or an electrical field. 

As long as one is taking the limit e ^ with t, s, L 
fixed, the polymer is correctly described by the linearized 
version of Eq. (|19|l . 

dtr^ = -r"l' + (/rl)' (22a) 
/" = -0'|, (22b) 

From the magnitude of the noise-correlations, given by 
Eq. I|10b|) . it may be inferred that ^± to leading order 
obeys 

(^±(s, mA-^\ t')) - 2{l/£j,)5{s - s')5{t ~ t') , (23) 

where we used kBTC,± ~* ^ choice of units and 

ly = Sij is the identity matrix. 

From Eq. (|22b(l . it is seen that the tension profile of 
a polymer that is forced at the ends is linear; slope and 
offset are fixed by the boundary conditions. The (higher 
order) longitudinal displacements are slaved to the trans- 
verse ones by the arc length constraint, Eq. H15() . At the 
present level of approximation the exact equations of mo- 
tion have reduced to a linear equation for the transverse 
displacements alone. 



A. Generalized transverse response 

In many practical cases, for instance if the poly- 
mer is symmetrically pulled apart by a (possibly time- 
dependent) force, the tension / = f{t) is to lowest order 
spatially homogeneous such that the equation of motion, 
Eq. I|22a|) . reduces to 

dtr^ = -rT + .fr'i+^± . (24) 

Let us anticipate at this point that we will identify an 
inherent length scale separation, Eq. (|58|l in Sec. be- 
low, according to which the tension can be considered 
as slowly varying in space, such that Eq. H24|) describes 
the polymer dynamics locally (and even globally at late 
times). As a consequence, the solution of Eq. H24|) for 
a given spatially homogeneous tension history f(t) be- 
comes an important ingredient of the nonlinear theory 
and shall be analyzed in the following. 

Formally, the linear Langevin equation Eq. (|25|l is 
solved by 

r±{s,t)= ds' dt' x{s,s'-t,t')i^{s',t') . (25) 

Jo J ~oo 

The Green's function x{s, s';t,t') satisfies 

dtx{s,s';t,t') = -dtxis,s';t,t') + fit)d^,xis,s';t,t') 
+5{s- s')5{t-t') (26) 

and appropriate boundary conditions. It may be inter- 
preted as a causal response function that describes the 
spreading and the decay of contour undulations induced 
by a transverse force impulse at location s' and (elapsed) 
time t'. Eq. H25|l therefore can be said to represent the 
conformation at time t in terms of the accumulated re- 
sponse to the transverse noise history ^±{s' ,t') along the 
contour. 

In the general case of a time-dependent tension, it 
can be quite difficult to determine the Green's function 
X that obeys the prescribed boundary conditions, be- 
cause eigenmodes and eigenvalues of the linear operator 
(/(t) -I- d^)dg depend on the value of f{t) and thus be- 
come time-dependent. In terms of Fourier modes, on the 
other hand, a translationally invariant Green's function 
x{s — s';t,t') can easily be found, below. As will be 
detailed in Sec. IIVI this function describes the universal 
bulk dynamics far away from the ends, while a correction 
term that manifestly breaks translational invariance has 
to be added "close" to the ends to correct for the actual 
boundary effects. 

To formalize this decomposition into bulk and bound- 
aries, the full response function x m^iy be written as a 
superposition 

X(s, s'; t, t') - x{s - s'; t, t') + xtc{s, s'; t, t') , (27) 

where x{s — s'; t, t') and Xfcc(s, s'; i, t') represent a trans- 
lationally invariant part and the boundary correction, re- 
spectively. The former is taken to satisfy 

dtx{s:t,t') = [-dl + f{t)]dlx{s-t,t') + 5{s)5{t^t') (28) 



6 



on an infinite arc length interval. With help of Fourier 
modes, 



/oo 
dsx{s;t,t')e 
-OO 



-iqs 



Eq. (EH|) reads 

dtx{q\ t, t') + \{q, t)x{q; t, t') = 5{t ~ t') 
where \{q^t) is the dispersion relation |55j| 



(29) 



(30) 



(31) 



By the method of integrated factors, the solution to 
Eq. |(2n|) is found to be [s^ 



x(g;i,<') = 0(i-Oexp 



- / dt\{q,t) 



(32) 



which may be checked by direct substitution. The real 
space susceptibility is given by the inverse Fourier trans- 
form of Eq. 



"dq 



X{q;t,t')cos{qs) 



(33) 



where it has been used that xily t') is even in q. 

The part Xtc(s, s'; t, t') of the susceptibility, which is 
not translationally invariant, satisfies the homogeneous 
differential equation 



dtXbc 



(34) 



and has boundary conditions that have to compensate 
for the generally inappropriate behavior of x(s — s'; i, <') 
at the boundaries. The difficulty of solving Eq. (|26|l has 
been shifted to Xbc{s, s';t,t'). However, for the time- 
dependent quantities to be studied below, this bound- 
ary term represents a relevant contribution only within 
a characteristic length £±{t) (defined below) close to the 
boundaries. For times small enough, such that £±{t) <C 
L, one may use Xbc derived on a semi-infinite polymer 
to approximate the situation near one boundary, say the 
one at s = 0. For simplicity, we will mostly refer to the 
model boundary conditions of "hinged" (h) or "clamped" 
(c) ends 01 j for which the full susceptibility on a semi- 
infinite arclength interval s G (0,cx)) is given by a sym- 
metric and antisymmetric combination of the bulk sus- 
ceptibility X, 

X| {s, s'; t, t') = x{s ~s';t, t') T x{s + s'- 1, t') . (35) 

Evidently, Xb. satisfies Eq. H26(l on s e (0,00), as well as 
the boundary conditions 

X/.(0, s'- 1, t') = = d^.XhiO, s'; t, t') (hinged) 
5,Xc(0, s'; i, t')^ Q = 9^x^(0, s'; i, t') (clamped) . 

A hinged end has vanishing transverse displacement and 
is torque free (vanishing second derivative), whereas a 
( "gliding" ) clamped end has a vanishing slope and is force 
free (vanishing third derivative). 



B. Local response of the stored excess length 

As a basis for our subsequent systematic analysis of 
tension dynamics and as an application of the above re- 
sults, we wish to determine the longitudinal motion im- 
plied by the linearized transverse stochastic dynamics, 
Eq. (|22a|) . To this end, consider a stiff polymer equi- 
librated under a constant tension /< at time zero, on 
which a spatially constant tension f[t) is imposed that 
varies deterministically for t > 0. We ask, how the den- 
sity of stored excess length g{s,t), defined in Eq. ifTI^ . 
changes in time by considering the ensemble average of 
the increase (during the time interval t) of the stored 
length 

Agis,t) = gis,t)~g{s,0). (36) 

This is an important observable, since it governs the lead- 
ing order contribution to the change 

Ai?||(i)^i?||(t)-i?||(0) 

in the projected end-to-end distance R\\(t), defined in 
Eq. (|16|) . which is can be directly measured in dynamic 
single polymer experiments p^ . Here, the (average) end- 
to-end axis of the polymer is assumed to be controlled 
by external means; e.g. by an external force field, flow 
field, boundary conditions, etc. To the relevant order, 
the precise measures taken to orient the polymer (strictly 
or on average) do not matter, and we find 



(Ai?||)(t) 



ds {Ag){s,t) + o{e) . (37) 



Here, the notation with the arguments s and t outside the 
brackets of {A.g) was introduced to emphasize that, even 
after averaging, these dependencies generally persist. 

For an explicit calculation of the stored length, we in- 
sert Eq. (|25|l into Eq. (|18|l and perform an ensemble av- 
erage upon employing Eq. (|23() 

/t pL pt pL 

dt' / ds' / di' / ds' dsx{s,s';t,t') 
-00 Jo J~oo Jo 

dsx{s,s';t,i'){Uis\t')-u{s',i'))/2 

—- / ds'dsx{-%s';t,t'f . (38) 

oo'-pV' ) Jo 



= 2 



For later convenience, we have allowed for a time- 
dependent persistence length £p{t) and an optional pre- 
stress /< ^ L^. The latter is also technically advan- 
tageous, since it acts as a physical rcgularization to sup- 
press modes with wavelength larger than the total length. 
It enables us to take the total length to infinity and 
to discuss the stored length (g) (s,t) on a semi-infinite 
arclength interval. For our ultimate goal of calculating 
(Ag) (s, t) an intrinsic rcgularization renders modes with 
wave length beyond a characteristic length scale £±{t) 
irrelevant, so that /< can eventually be set to zero if 
required. 
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Inserting Eq. H35() . valid for hinged/clamped boundary 
conditions, into Eq. (|38|l (with L ^ oo) yields 



Ag.) = f ^[g{q,t)-g{q,0)][l±cosi2qs)] 



g{q,t) = 2q' / dtx\q;t,t)/ep{t) 



(39) 



The general expression Eq. (|39|1 is now specialized to the 
scenario 



tension persistence length 

^ < : /< = const. (p = const. (40) 

t>0: f{t) ip/e = const. 



As preparation for a more general discussion, we have 
included the possibility of a sudden change in persistence 
length by a factor l/6'>0att = 0. 

With a constant tension at negative times, the time- 
integral in Eq. H39|l can be evaluated from i ~ — oo to 
t = 0, after which we obtain 

Sil't)£p = ^p^ + 2dq' fdtt{q:t,i). (41) 



+ f< JO 

Therewith, the full expression for the average change in 
stored excess length density becomes 



AgH)is,t) = / 



1 



Trip + /< 
Jo 

[l±cos(2gs)] , 



1 



(42) 



where F{t) is the integral of the tension over positive 
times, 



Fit) = / difii) 



(43) 



One observes that the integral in Eq. (|42|l is well de- 
fined, even for /< ~ 0, because the integrand vanishes for 
g — > 0, oo. It is dominated by wave numbers for which 
the exponents become of order one, i.e., the dominant 
modes have wave numbers q for which the characteristic 
relaxation time 



T, = {q' + q^F/t)-' 



(44) 



is of the order oft. This suggests to define a characteristic 
length £±{t) as the wave length for which the relaxation 
time is just t, 



l^£l^ilh + F{t)] . 
Asymptotically £_i_(t) is given by 

t^/" , for i < (F/t)-^ 



£^{t) 



Fi/2 , for i > (F/t)-^ 



(45) 



(46) 



TABLE I: The transverse equilibration length £±{t) and the 
tension propagation length £^^ (t) both exhibit a crossover at a 
time tf = f which depends on the external force f (here, 
for the Pulling problem with f ^ L^'^)- 





£±{t) 


£\\{t) 


t < tf 


^1/4 


t'^'ii^/cy^^ 




t > tf 


^1/2^/2 


t'/*f'/'(^p/C)'/' 


]£ 



Due to the competition between bending forces (cx 
'''-l£±^) and tension (cx rj_{F/t)£^^), the growth of £± 
thus exhibits a dynamic crossover from free relaxation 
{£j_ ~ t^/**) to relaxation under tension {£j_ ~ ^/F) at a 
characteristic time t-j = {F/t)~^ (Tabled/left for a con- 
stant tension equal to the external force f ). 

As we will explicitely demonstrate in the next section, 

the change (^Agn^ in stored length saturates at a con- 
stant value for distances to the boundaries much larger 
than the characteristic length £±{t). This "bulk" value 
is given by js^l 

iA-g)it) . r^|^(e-'^''^--^*)l-l) 
Jo '^£p 1 9 + /< ^ ^ 

+ 29q' /We"2«'[''(*-*>^(*)-^(*)]l .(47) 



The quantity (Ap) (t) will be central in our systematic 
analysis of tension propagation and relaxation, because it 
turns out to determine the local curvature of the tension 
profile. Each of the two terms inside the curly brackets 
of Eq. 1)47(1 have a direct physical interpretation. Since 
the parameter 9 tunes the strength of the thermal kicks, 
it is seen that the ^-independent first term represents the 
deterministic change in the excess length that is stored 
in mode q in absence of any stochastic force. For pulling 
forces _F > its sign is always negative, since both the 
internal elastic and the external driving forces act to 
straighten the filament. On the contrary, thermal kicks 
represented by the (strictly positive) second term are fa- 
voring undulations. 



IV. BREAKDOWN OF ORDINARY 
PERTURBATION THEORY 

The previous sections employed "ordinary" perturba- 
tion theory (OPT) in the small parameter e, leading to 
a linear equation of motion to lowest order. As detailed 
below, the use of OPT is, however, limited to long times 
even for e ^ 1. The predictions derived above for the 
longitudinal segment motion turn out to be incompati- 
ble with the longitudinal force balance on short times. In 
particular, Eq. (|47|l reveals an infinite longitudinal fric- 
tion for t 0. This section extends a heuristic argument 
of Ref. to resolve this problem. The following, further- 
more, elucidates a very general feature of the non-linear 
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response, namely, a crossover from "weak-" to "strong- 
force" behavior. Finally, it reveals the crucial length- 
scale separation underlying our subsequent systematic 
analysis. 

The breakdown of OPT becomes evident when we try 
to use Eq. (|42|l to evaluate the longitudinal segment mo- 
tion in a non-equilibrium situation. For dcfinitcness and 
as a telling example, let us consider an initially equili- 
brated polymer that is suddenly pulled longitudinally by 
a constant force f at both ends, i.e., at time t = 0, the 
tension at the ends is suddenly increased from to a 
given positive value f , 

f = /(O, i > 0) = f(L, t > 0) {Pulling-scenaiio). 



As a consequence of Eq. (j22b|l . the tension / = f is 
to lowest order spatially uniform and fixed by the driving 
force f at the boundaries. Our above result for the change 
in stored length due to a given tension history, Eq. (|42|l . 
thus applies and evaluates to 



dq 



f 



nl^q^q^ + j) 
X [1 ± cos(2gs)] . 



After the variable substitutions k = qf (J 
and T = \^t = t/tf , Eq. gHIl takes the form 



(48) 
(49) 



with the scaling function 



Eh (cr, t) 



1 



TT fc2(fc2 + l) 

X [1 ± cos(2A:cr)] . 



- 1 



(50) 



Another variable substitution k = kr^^^ generates factors 
k^ + k^T^^^ in exponent and denominator, which can be 
replaced by k^ in the asymptotic limit t 1. Just as 
for the dispersion relation Eq. H44() . tensile forces oc k^ 
may be neglected as compared to the dominant bending 
forces cx fc"* for times smaller than the crossover time if . 
In the opposite limit r ^ 1, the reverse approximation 
applies. After a variable substitution k = kr^^^ , factors 
k^T^^ + fc^ appear, which may be replaced by fc^. We 
thus find that the two-parameter scaling form Eq. H5U|) 
collapses onto one-parameter scaling forms for small and 
large times. 



E. ~ -r3/4s^(aT-i/4) 



-1/2 





oo 



with scaling functions given by 

'^dk 1 - e-2fe' 



s|(C) 

mo = I 5 



dk 1 - "-^fe^ 



k^ 



1^ A;2 



l±cos(2fcO , (51) 
l±cos(2A;^)j . (52) 




FIG. 3: The absolute value 



AE1 



of the spatially varying 



part AE'J = E'J; (^) — E'J (oo) of the scaling functions E'^; 

defined in Eqs. 1511 152II . which happens to be the same for 
hinged/clamped boundary conditions. 



Note that the spatial part of these scaling functions de- 
picted in Fig. O decays to zero within several rescaled 
time units. As an important consequence, we note that 

the part of {^Agh^ that depends on the boundary con- 
ditions really only matters close to the boundaries, i.e., 
up to a distance for which the scaling variable becomes 
of order one. In fact, this distance can be identified with 
(<) , as defined in Eq. (|45() , because the scaling variable 
of E is given by crr-^/* = s/t^''^ = s/e±{t) for r < 1 
and ar^^^^ = s/i^{t) for t ^ 1, respectively. 

The bulk of the polymer stores length according to 
the universal part (Ag) [t) (independent of the boundary 
condition), which asymptotically takes the form 



(A^) (0 



r(7/4)-if €;i(t/2)3/4 m , for t < if 
i(2ft/7r)i/2 , fort>tf 



(53) 

For alternative derivations of the short-time linear re- 
sponse law ^ t^/'', see Refs. |^ 0, The crossover 
time = f^^ is the time where the external force f 
equals the Euler buckling force i']^{t) corresponding to 
the correlation length it) . 

These results comprise the predictions of ordinary per- 
turbation theory (OPT) to leading order. As evident 
from Eq. H24() . longitudinal friction forces have thereby 
been completely neglected, because they are of higher 
order in e. However, on the semi-infinite arclength inter- 
val considered here, a spatially constant change in stored 
length - no matter how small - implies via Eq. H37() an 
infinitely fast change dt (Ai?||) of the longitudinal exten- 
sion and thus friction. Hence, OPT must fail on suffi- 
ciently large length scales. Let us define i*(i) as the 
length scale beyond which OPT breaks down for a given 
time t. In the present case of a suddenly applied pulling 
force, this critical length can be estimated from the phys- 
ical requirement that the total longitudinal friction is not 
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only finite but at most equal to the driving force f , i.e., 

f ~ (Ldt (Ai?||) ^ -QL / dsdt {Ag) {s,t) . (54) 
Jo 

For an average stored length given by the spatially con- 
stant value Eq. H53() , this condition is met if the polymer 
length is smaller than 



for t <s:tf 



Vc(^^)^/^ fort»tf 



(55) 



OPT can thus only be valid on length scales (much) 
smaller than Li,{t). The condition 



(56) 



implicitly defines the critical time scale above which 
OPT applies to the whole polymer, and below which 
OPT is limited to subsections shorter than i*(t). 

In summary, the omission of tension propagation has 
been identified as the reason for the breakdown of the 
OPT predictions. On a heuristic level |^, the problem 
can thus be resolved by requiring that the applied ten- 
sion is not immediately perceptible everywhere in the fil- 
ament, but propagates a finite distance i\\{t) during time 
t from the ends towards the bulk of the filament. Hence, 
only segments up to a distance £|| it) from the ends are set 
into longitudinal motion. If the length t\\{t) over which 
the tension varies is smaller than the critical length Li,{t), 
it is ensured that the longitudinal friction does not ex- 
ceed the driving force. In heuristic discussions, it was 
generally assumed that both lengths can be identified up 
to numerical factors of order one, 

^11 (i) ~ Li,{t) , t -^ti, , (heuristic hypothesis) (57) 

as summarized in Table ^^right. The scaling assump- 
tion in Eq. H57|) turns out to reproduce the correct scal- 
ing for most cases (with, however, interesting exceptions 
elaborated in Part II ^ "^^11 as, in Ref. 0|). The 
"weak- and strong- force" limits l\\ (x t^^^ and £|| oc t^^^ 
of Refs. 13,0 are with Eqs. (|551[S7|I recovered as "short- 
and long-time" asymptotics. The crossover at time tj sig- 
nals the change from "free" to "forced" relaxation and is 
inherited from the one of the scaling function S, defined 
in Eq. (EOJ. 

Finally, we note that the tension may be considered as 
"slowly" varying in space because £|| (t) oc 0{e~^/^) is for 
e ^ larger than any length that does not dependent 
on the small parameter (for i, f, L fixed). As it turns 
out, the most important length of the latter type is the 
dynamic correlation length £±_ (t) for transverse displace- 
ments. Namely, the scale separation 



indicates that the tension is nearly constant on the equi- 
libration length scale i±{t) for transverse displacements. 
Intuitively, it should be clear that this simplifies the fur- 
ther analysis considerably, because it allows to apply 
(certain) results locally that are derived for spatially con- 
stant tension. Formally, Eq. (|58|l lends itself as a starting 
point for a multiple- scale calculus, which separates the 
physics on different scales to obtain an improved pertur- 
bation expansion that is regular in the limit t — > while 
e <C 1 is fixed [H^- The procedure, detailed in the next 
section, is similar in spirit to the procedure for athermal 
rod dynamics [4^ , but some complications related to the 
stochastic nature of the equations of motion have to be 
faced. The final result will be an effective deterministic 
description of the tension on the macro-scale i\\ (t), where 
the stochasticity on the micro-scale £±{t) has been inte- 
grated out. 



V. MULTIPLE SCALE ANALYSIS 

We introduce a rapidly and a slowly varying arc length 
coordinate, x = s and y = se^, respectively, where the ex- 
ponent 7 > will be fixed later. The dynamic functions 
r± and / are now considered to depend on both vari- 
ables {f,r±} {f{x,y),r±{x,y)}, where x and y are 
treated as independent. The original arc length deriva- 
tive of those functions then becomes 



ds = dx\y + e'^dyla; . 



(59) 



The noise ^ = 0{e^^^) being the source of any trans- 
verse displacements suggests an expansions of the dy- 
namic variables r± and / in powers of e^^^, 

r^{x,y) - ei/2/ii(a;,y)+o(ei/2) , 

f{x,y) = fo{x,y) + e^/'^fi{x,y) + ef2ix,y) 

+o{e) . (60) 

In the case of isotropic friction (i.e. C_l — C||)> the stochas- 
tic forces have no intrinsic scale. Hence, they can only 
depend on the microscopic variable, ^ = e^/^^i(x). In 
the anisotropic case, the friction forces and, hence, the 
stochastic forces, are coupled to the orientation of the 
filament, so that one has to assume a power expansion 



|^(a;,y) = £i/2|^,i(.T)+o(ei/2 

(x, y) = e^'^i\\s{x) + e^||,2(a^, v) + o(e) 



(61) 



(58) 



The y-arguments in Eq. (|61|l are inherited from the 
2/-arguments of r' entering the noise correlations in 
Eq. |T0B|) via the friction matrix, defined in Eq. (jH)). 
The y-dependence would disappear for isotropic friction. 
Note, that the leading order ^i(a;) still depends on the mi- 
croscopic variable x only, because the anisotropy merely 
enters the higher orders (a formal argument is given in 
Add.IbTIi. 
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In the following, it is crucial to require that the ex- 
pansion coefficients in each order are bounded, so that 
we obtain a uniformly valid power expansion |48l | in 
terms of the small parameter e. Inserting all expansions, 
Eqs. I5n|l . into the equations of motion, Eq. (|19al 
ITflfe)! . yields 

= [d,h, + d% - d^od.hi) - l^^i] 

+ o(ei/2) (62a) 

= dlfo + e^2d,dyf„ + [d^h + 

+ e-^+^'^dyd^h + o(e; e'^) ■ (62b) 

In order to arrive at Eq. (|62bp we used the local arc length 
constraint, Eq. \ib\ . By 

c. 



^2(x,y) = ^dt{d.,hif + [l - C)dA{d.hi){dthi)] 



\dt{d,h,f -\dl 



fo {dxhi 



(63) 



we have summarized terms nonlinear in hi. The first 
term in X2, which is proportional to C, accounts for the 
longitudinal friction and is thus responsible for the short- 
time divergence encountered in the heuristic discussion of 
Sec.HVI 

The 0(1) part of Eq. l|62H)l . d^fn = 0, together with 
the requirement of fo being bounded for large x, implies 
that 



/o(a;,y) = /o(y) 



(64) 



is independent of x. Hence both the 0(1) and the 0{e'') 
term of Eq. (j62b|) vanish. Requiring the 0{e^^'^) coef- 
ficient to be zero fixes /i up to an integration constant. 
The value of /i does not affect the evolution of hi because 
it does not enter the 0{e) coefficient in Eq. (|62a|l . Thus, 
the precise value of /i does not change the pathologi- 
cal behavior of longitudinal friction oc (^dt{dxhi)^ , which 
is why we shift the discussion of fi to App. IB 21 From 
the latter we merely need that dydxfi{x,y) = which 



renders the 0{e'<+^/'^)-term in Eq. l|62B|l zero. Then the 
next higher order is either 0{e^^) or 0{e) depending on 
the value of 7. With Eq. we can solve the (e^/^) 
part of Eq. (|62a|) for hi {x, y) along the lines of Sec. IIII Al 
and use the result to evaluate X2{x, y). It then turns out 
that the first term in X2 (the longitudinal friction) would 
require /2 to grow without bound with increasing system 
size to render the C'(e)-expansion coefficient in Eq. Ij62b|) 
finite, if they were required to balance each other. This 
represents the same unphysical divergence that is respon- 
sible for the breakdown of ordinary perturbation theory, 
discussed in Sec. IIVI 

In order to obtain an improved perturbation theory, 
we attempt to balance the nonlinear term by the 0{€^^) 
term after choosing 7 = 1/2; i.e., the exponent 7 is fixed 



such that the expansion coefficient /2 remains finite in the 
semi-infinite system considered here |59j . The equation 
fixing /2 thus reads 



(65) 



Given hi and /o the last equation can be solved for /2 



X nX 



h{x,y) = j dx j dx \-dlfo{y) 



-X 



2(5,2;) -C|l,2(^)} 



(66) 



For /2 to be bounded for large system sizes, we have to 
require 



dlUv) = -X2-il, iy) 



(67) 



where the over-line denotes the spatial average over the 
rapidly varying coordinate x, 



Tiy) 



lim / ^9{x,y) 



(68) 



for a function g{x,y). The expansion coefficient /2 would 
show a divergence quadratic in the system size if Eq. H67() 
was not satisfied. Hence, the y-dcpendcncc of /o(y) must 
be fixed such that the expansion coefficient /2 remains 
finite. For a finite polymer, the limit Z 00 is not to be 
taken literally, though. Rather, the average in Eq. H68|) 
is required to become independent of I to leading order 
in e for / much smaller than the system size L. 

As it turns out, the only quantity in Eq. (|67|1 that 
does not vanish upon .x-averaging is the first term in 
X2, the longitudinal friction. This is easily seen for 
all other terms in X2 and the /j-tcrm. They are total 
derivatives with respect to x of products of expansion 
coefficients that are (required to remain) bounded (non- 
secular by definition. Hence, the x- integrals of those 
total derivatives are bounded and the x-averages vanish 
upon formally taking the coarse-graining length / — > cx) in 
Eq. (|68() . The noise term also represents a total derivative 
with respect to x. The average of that term represents a 
stochastic variable with an amplitude that scales as 1/Z 
and, hence, also vanishes in the limit / — *■ cxd. 

Dropping all terms that vanish under coarse-graining 
and integrating over time, Eq. H67|) takes the form 



-^dlFoiy) 



[dxhif {y,Q)-{dxhiY {y,t) 



(69) 



where Aq{s, t) is the change in stored length, as defined 
in Eq. H36|l . of a semi-infinite polymer for the tension 
history 



F{t)=Fo{y,t) 



(70) 



Note that the dependence on the slowly varying arclength 
coordinate y enters the tension history in Eq. H70|) only 
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paramctrically. The same holds true for the calculation 
of the right hand side of Eq. (|^ . 

The closed set of equations Eqs. (1691 162a|l represents 
the lowest order of the multiple scale perturbation ex- 
pansion. It incorporates the feedback mechanism already 
found in the heuristic discussion above. The evolution of 
transverse displacements implies longitudinal motion via 
the arc length constraint, and, according to Eq. (|69|l . the 
corresponding longitudinal friction sets the polymer un- 
der tension. This, in turn, feeds back onto the evolution 
of transverse displacements, Eq. (|62a|l . typically acting 
as to reduce the longitudinal friction. 

For solving Eqs. (|69l I62a() self-consistently, it would 
be handy to perform an ensemble average (...) on the 
right-hand-side of Eq. (|69|l . because we could then apply 
expression Eq. H42() for (Ag). We argue that such an en- 
semble average is indeed justified, because Ag {t) is, in 
fact, a deterministic quantity as a consequence of the cen- 
tral limit theorem. Recall that the quantity {Ag) (s, t) is 
dominated by transverse modes of the wave length £± (t) 
given by Eq. (|46|l . Since wavelength much larger than 
ij_(t) are not relevant in the mode sum, the dynamic 
length i± (t) can be interpreted as the correlation length 
of the change in stored length, i.e., the length over which 
the correlation function Cs,t(z) = {Ag{s,t)Ag{s + z,t)) 

varies. Hence, the integral Xi = J^^ds Ag{s,t) may be 
understood as the sum of l/£±{t) weakly correlated ran- 
dom variables. For large I, the distribution function of Xi 
thus becomes Gaussian, and the variance of Xi grows lin- 
early with the number of independent random variables, 
({Xi — (Xj))^) cx L As a consequence, the distribution 

of the average Xi/l — > Ag'' approaches a delta function 
as Z — > oo. 

An "additional" ensemble average therefore does not 
change the value of the right hand side of Eq. (|^ . 
Evaluating the spatial average after the ensemble aver- 
age levels out the boundary term of expression Eq. 

for (^Agh^ and reduces it to its bulk value (Ag) (t), de- 
fined in Eq. (|47|l . which is completely independent of the 
boundary conditions. We thus have the important rela- 
tion 

A^\t) = (^Agj,)\t) ^ (Ag) (t) . (71) 

We would like to emphasize, that our argument for re- 
placing the spatial by an ensemble average requires a 
finite driving force, such that limj^cxD Ag approaches a 
finite value. This specifically excludes the linear response 
limit, i.e., the limit of vanishing external force f ^ 
while e ^ 1 is fixed. In this case, the tension dynam- 
ics has to be described by Eq. Ht)9|) . which is stochastic 
even to leading order (we will come back to this point in 
Part II). 

Given the external driving is finite such that Eq. 1171(1 
may be applied, Eq. (|69|l takes the form 

92f(s, t) = (A^) [Fis, t < i), t] , (72) 



As 




s — r 



FIG. 4: Tensile forces acting on a polymer subsection of size 
As. Balancing these forces with the drag that arises from 
the longitudinal velocity of magnitude dtr^^{s,t), one esti- 
mates AsC||9tr||(s,t) f» /(s - As/2) - /(s + As/2) in the 
weakly bending limit. For infinitesimal As, this becomes 
C||c'tr|| (s, t) « —dsf{s,t). A further spatial derivative yields 
C,\ldtdsr\\{s,t) ~ —dsf{s,t), which is the time derivative of 
Eq. 1721 (in original units) up to a spatial and ensemble aver- 
age on the left-hand side, which correspond to an adiabatic, 
respectively, equilibrium approximation. 

where we introduced s ~ ye^^^ ~ x again and made the 
parametric dependence of (A'g) on the tension history ex- 
plicit. The deterministic tension dynamics, as described 
by Eq. (|7^ , provides the sought-after rigorous local gen- 
eralization of the heuristic argumentation of Sec. IIII Bl 
local longitudinal motion is driven by tension gradients 
(like in a thread pulled through a viscous medium). 

Upon inserting our result in Eq. H47(l for the right- 
hand-side, we obtain a nonlinear, partial integro- 
differcntial equation (PIDE) for the tension history 
Fis,t), 

We have arrived at a closed description of the polymer 
dynamics to lowest order in MSPT that consists of two 
parts. On a length scale t±{t), Brownian motion gives 
rise to fiuctuations of transverse displacements that are 
described by the linear Eq. (|62a|l . This stochastic dif- 
ferential equation, in turn, adiabatically depends on a 
tension profile that varies on a much larger scale i\\ (t) 
and satisfies a deterministic nonlinear equation of mo- 
tion, Eq. lO. 



VI. CONCLUSION 

We would like to conclude the present general discus- 
sion of tension dynamics with a simple physical interpre- 
tation of the outcome. Eq. (|72II73|I . of our multiple scale 
analysis. Effectively, our MSPT analysis is a rigorous 
justification of certain approximations that can be made 
to analyze the tension dynamics in the small e limit. Ac- 
cording to Eq. (|7^ . the curvature of the integrated ten- 
sion is (up to a constant) given by the ensemble average 
of the local stored length release. As shown in Fig. 01 
it may be conceived as (1) a force balance equation be- 
tween the locally acting tensile and longitudinal friction 
force, in which (2) the latter may be computed from the 
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equations of motion for the bulk of an equilibrated poly- 
mer under a spatially constant, though time-dependent, 
tension. The longitudinal friction force is obtained via 
taking a time-derivative of this coarse-grained dynamical 
force extension relation, Eq. l(17|l . The first order MSPT 
equation of motion neglects longitudinally acting bending 
forces, employs an adiabatic approximation and assumes 
local equilibrium. The latter fully retains memory effects 
and therefore must not be mistaken as an approximation 
of quasi-stationary dynamics, which, in Part II [To| . will 
turn out to be a valid approximation only for specific 
driving forces in a particular time regime. 
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APPENDIX A: ELASTIC FORCES 

From the effective free energy of a worm-like chain, 
Eq. ||2J), we seek to determine the elastic force goi(s) per 
arc length that a polymer of a given conformation r(s) 
exerts at arc length s onto its surroundings. To this end. 



TABLE II: Some important notations 



Symbol(s) 


General Meaning 


L 


total length of the worm-like chain 




K 


bending stiffness 




e 


small parameter, defined such that = 




r±{s,t) 


transverse displacement; r±_ = 0{e^^^) 




r\\{s,t) 


longitudinal displacement; = 0(e) 




g{s,t) 


stored- length density; g = r'l/2 + 0(6^) 


= 0(e) 


R 


end-to-end distance 






end-to-end vector projected onto the Ion 


J. axis 




equal up to numerical factors of order 1 




(X 


proportional to 
asymptotically equal 




e 


persistence length 

Ip changes by a factor 1/S>0atf = 0. 






thermal energy 




fis,t) 


line tension 






equilibration scale for transverse bending modes 




scale of tension variations 




f 


external force 







external force per arc length 






thermal force per arc length 





let us assume the polymer was subject to a constant ex- 
ternal force field equal to —gei{s), the negative of the 
local elastic forces. Then, the conformation f(s) is in 
balance with the external force and minimizes the total 
free energy 



(Al) 



with 



T^extH = / dsg,i{s)r{s) , (A2) 

JQ 

under all possible paths that obey the local inextensi- 
bility. Eq. Hence, requiring a vanishing free energy 
change 6Ti. for an infinitesimal (permitted) change Sr in 
the space curve leads to the sought-after relation between 
the elastic forces and the contour. The central question 
of the minimization problem is how to deal with the local 
constraint. 

Here, we show that the common |49j introduction of a 
Lagrange multiplier function ensuring the local inextensi- 
bility constraint is not necessary. We will instead present 
a minimization procedure that considers only variations 
that obey the incxtensibility constraint to leading order. 

Consider a test contour 



r(s) = r(.s) + Sr{s) 



(A3) 



which is infinitesimally displaced by dr{s) from the equi- 
librium contour r{s). The inextensibility constraint, 
Eq. (^3), is fulfilled to 0{dr) if we only consider displace- 
ments that are constructed from another infinitesimal 
vector field e(s) according to 



Sr'{s) = e{s) X r'(s) 



(A4) 



Those displacements are transverse to the local tangent 
vector, so that (5(r'^) is quadratic in Sr. They correspond 
to local rotations of the tangents. 

The variation of the contour induces a variation STi. of 
the total elastic energy H = Hwlc + Wcxt of the form 



5n 



[ ds {Kf"5r" 
Jq 



- [ ds 
Jo 



^9ciSr) (A5a) 
f J dsgci^ Sr' + b.t.(A5b) 
- /oi ) (e X r') + b.t. (A5c) 



- [ ds [nf" 

- [ dse[f' X (Kf"'+/el)]+b.t. 

Jo 



(A5d) 



Here, we performed a partial integration to obtain 
Eq. ljA5bp . which introduces some boundary terms ab- 
breviated by "b.t.". In the subsequent step we inserted 
Eq. HA4|) for Sr' and introduced the arc length dependent 
force 



/oi (s) = / dsgoiis) 
Jo 



(A6) 
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Finally, we used the property A ■ {B y. C) = C ■ {A y. B) 
of the triple scalar product to obtain Eq. ljA5d|) . 

If r(s) is indeed the equilibrium contour then STi, has 
to vanish for all variations parameterized by e{s) and 
Sr{L). Therefore, the term in the square brackets of the 
integrand in Eq. ljA5d() has to vanish, 

r' x{Kf"' +f,i)=0 , (A7) 

so that we recover Eq. Q , used in Sec. ^ In addition, 
the boundary terms 

b.t. = K{r"dr')\^ +6r{L) ^dsg,, 

Jo 

= Mr"(exf')]|o^ + [/ei<5r]|^. 

have to cancel, implying the requirements 

/ci(i) = (A9) 
ir' X f") \o,L = . (AlO) 

The condition expressed by Eq. (|A9|) simply states that a 
force balance can only exist if the external forces sum up 
to zero. Using the inextensibility constraint, Eq. (Q, it is 
seen that f" • f' = 0, so that Eq. (|A10p can be rewritten 
as the boundary condition 

r"\ox = . (All) 

The curvature, which is proportional to the local torque, 
has to vanish at the (free) ends. 

Note that the over-bar in f(s) to denote the equilibrium 
contour is dropped in the main text, for simplicity. 

APPENDIX B: MULTIPLE SCALE ANALYSIS 
(DETAILS) 

1. Thermal forces 

In the multiple scale perturbation theory, presented in 
Sec. one should, in principle, assume an expansion 
^{x,y;t) = e^^^^i{x,y,t) + o(e) where the leading order 



noise is a function of the coarse-grained variable y. 
However, from the fundamental correlations obeyed by 
the original noise function, Eq. HlObp . we have to require 
to leading order 

{Us,es;t)^i{s' ,es';t)) = 2(I/L)S{s - s')S{t - t') (Bl) 

for all e = L/£p. Apparently, the right hand side of 
Eq. (|ljl|l does not depend on e. As a consequence, the 
left hand side cannot depend on e either, in particular 
not on es' or es. Hence, the two-point correlations are 
independent of the slowly varying arc length coordinate. 
Using Wick's theorem, we may argue in the same way 
for all higher-order correlation functions as well and con- 
clude that the leading order stochastic force is itself 
independent of y. 

We note, that the inverse length appearing on the right 
hand side of Eq. (|B1|) is due to the definition e = L/£p 
of the small parameter for a stiff polymer. In the case of 
a strongly pre-stretched polymer, the small parameter is 
defined as e = (£pf )~^/^, see Eq. (fT^ and the subsequent 
paragraph, so that L^^ in Eq. (jBlj) has to be replaced 
byf-i/^ 

2. Next to leading order tension 

Requiring the 0{e^^^) coefficient in Eq. Ij62b|l yields a 
noisy first order tension 

fi{x,y) ^ / [^|i,i(0)-^||,i(i)] -f 6i(y)a;-f ai(?/). 
Jo 

(B2) 

For /i to be bounded in x, the term bi{y)x has to cancel 
the linearly growing the noise-term on the right hand 
side, 

bi= lim - Cdi [e||.i(0)-e||,i(5)] . (B3) 

However, important for the multiple scale analysis in 
Sec. is merely that bi is independent of y, so that 
dydxfi{x,y) = 0. 
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